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Abstract 



This paper describes the use of the time-of-flight (TOF) technique as a particle 
identification method for the HERMES experiment. The time-of-flight is measured 
by two 1 x 4m 2 scintillation hodoscopes that initially were designed for the first-level 
trigger only. However, the suitable time structure of the HERA electron beam allows 
an extension of their functions to also measure the TOF for low momentum hadron 
identification. Using only these conventional hodoscopes, good particle identification 
was achieved for protons and pions in the momentum range up to 2.9 GeV/c and 
for kaons up to 1.5 GeV/c. 

PACS number 29.30.-h, 29.40.-n, 29.40.Gx, 29.40.Mc 
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1 Introduction 



The HERMES experiment at DESY [1] is a 
second generation polarised deep-inelastic 
scattering experiment to study the spin 
structure of the nucleon. Several experi- 
ments over the last decade have provided 
accurate data on the polarisation asymme- 
try of the cross-section for inclusive scat- 
tering where only the scattered lepton is 
detected. Further knowledge of the origin 
of the nucleon's spin can be gained by study- 
ing semi-inclusive processes involving the 
detection of hadrons in coincidence with 
the scattered lepton. This increases the de- 
mands on hadron identification of the de- 
tection system. 

The time-of-flight (TOF) particle identifi- 
cation method is a fast, inexpensive and 
efficient technique [2] for hadron identifi- 
cation. Its implementation at HERMES is 
possible because of a) the presence of two 
scintillator counter walls HI and H2 in the 
HERMES spectrometer, and b) the fine 
time structure of the HERA electron beam 
(with bunch lengths of 27 ps and time be- 
tween bunches of 96 ns) , which allows par- 
ticles from different bunches to stay com- 
pletely separated in time. 



2 Counter Design 



HERMES is a forward spectrometer with a 
large dipole magnet, a set of tracking de- 
tectors, and particle identification (PID) 
detectors consisting of a Cerenkov detec- 
tor, a Transition Radiation Detector 
(TRD), a preshower detector and a calori- 
meter [1]. In 1998 the threshold-Cerenkov 
counter was upgraded to a Ring Imaging 
Cerenkov (RICH) detector [3]. A scintilla- 
tor hodoscope (HI) and a Pb-scintillator 
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Fig. 1. Schematic view of particle tracks and 
hodoscope paddles (upper part). 

preshower counter (H2) (Fig. 1) provide 
trigger signals and particle identification 
information. Both counters are composed 
of vertical scintillator modules (42 each in 
the upper and lower detectors), which are 1 
cm thick and 9.3 cm x 91 cm in area. The 
material for the modules is BC-412 from 
Bicron Co., a fast scintillator with large 
attenuation length (300-400 cm for scintil- 
lation light). The scintillation light is de- 
tected by 5.2 cm diameter Thorn EMI 9954 
photomultiplier tubes coupled via a light 
guide to the outer end of each scintillator 
(away from the beam plane). The mod- 
ules are staggered to provide maximum ef- 
ficiency with 2-3 mm of overlap between 
each unit. Each hodoscope photomultiplier 
(PMT) signal is passively split with one 
output going to a LeCroy 1881M ADC and 
the other going to a LeCroy 3420 Con- 
stant Fraction Discriminator (CFD). The 
individual CFD outputs are fed to LeCroy 
1875A time to digital converters (TDC), 
which measure the time-of-flight using the 
HERA-clock reference signal. This sig- 
nal corresponds to the moment when the 
HERA bunch cros-ses the center of the tar- 
get. When some interesting event configu- 
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ration is recognized by the HERMES trig- 
ger logic, the HERA-clock signal is enabled 
to start the TDC modules in common-start 
mode, with the STOP signals coming from 
each scintillator. The time base for the 
TDC is 50 ps/channel. 



the average speed of the electrons to the 
speed of light for each paddle individually. 
A single-counter time response is described 
by 

t = t pf + t lf + t (1) 



3 Calibration procedure 

The calibration procedure is based on the 
fact that electrons above 10 MeV are mov- 
ing at essentially the speed of light. Any 
measured deviations from this must be ar- 
tifacts of the experiment that should be 
corrected. 

In Fig. 2 the time distribution for elec- 
trons over the paddles is shown for part 
of one detector. The electron sample is se- 
lected by a combination of cuts on TRD, 
Cerenkov, Preshower and Calorimeter. The 
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Fig. 2. Non- equalized time dependence over 
X-coordinate 

time distribution is spread within each pad- 
dle as well as between paddles. The differ- 
ences between the paddles (due to different 
cable lengths, etc) were removed by setting 



where t p f is the particle time-of-flight from 
the interaction point to the scintillator 
plane, is the time needed for the light 
created in the scintillator to reach the PMT, 
t is the constant time offset specific to 
each paddle, arising from the PMT response 
time, cable delays, TDC calibration inter- 
cept, etc. The path length In of each track 
to hodoscope H was calculated using the 
reconstructed track parameters, along one 
straight line segment from the interaction 
point in the target to the mid-plane of the 
spectrometer magnet, and continuing along 
another straight line from there to the plane 
of the hodoscope H. This approximation is 
adequate because the magnet bend angle 
is less than n/20. The calibration for each 
paddle n of hodoscope H was done with 
electron tracks by fitting the y-distribution 
oil/v = t/ln shown in Fig. 3 to the known 
value of 1/c using a fourth-order correction 
polynomial F n (y): 
F n (y) = l/c-t/l H . 

The constant term in the polynomial in- 
corporates to in Eq. 1 for each paddle. A 
similar method has been used in [4] where 
a more detailed explanation for each coeffi- 
cient of the polynomial is given. Occasional 
shifts in the HERA-clock signal derived 
from the electron beam accelerator system, 
shown in Fig. 4, were compensated by re- 
fitting the constant term for each run of 
approximately 10 minutes. The final cor- 
rected distribution of l/v e — t/ljj + F n (y) 
for electrons is shown in Fig. 5. The reso- 
lution extracted from a Gaussian fit corre- 
sponds to <x=0.49 ns. 
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Fig. 3. l/v e non-equalized dependence on the 
Y-coordinate for electrons (single paddle). 
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Fig. 4. Y corrected l/v e dependence on run 
numbers (for electrons). 

4 Hadron Identification via 
Time-of-Flight 



From the relativistic momentum of the par- 
ticle 

p = m • /3/V(l-/3 2 ), (2) 
the mass m is extracted using the speed 
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Fig. 5. l/v e for electrons in hodoscope HI 
(equalized and run-corrected). 

f3 — v/c obtained through the hodoscope 
timing: 



m 



p 2 ( — 



(3) 



The squared mass of the particle was cho- 
sen as the parameter for the identification. 
To check the algorithm and calibration, 
the proton and kaon squared masses have 
been extracted giving values of m 2 = 0.88 
(GeV/c 2 ) 2 and m\ = 0.25 (GeV/c 2 ) 2 , re- 
spectively, which agrees well with the ex- 
pected values (Fig. 6). At momenta below 
2 GeV/c the kaon flux is two orders of mag- 
nitude smaller than the pion flux, so that 
no separation is possible via TOF for mo- 
menta above 1.5 GeV/c. Therefore in the 
momentum range 1.5 < p < 2.0 GeV/c, 
kaons are included in the pion spectrum 
as a negligible contamination {<!%). At 
higher momenta, a RICH detector takes 
over. 

Fig. 7 shows the squared hadron mass dis- 
tribution for the momentum region 1.5 < 
p < 2.0 GeV/c. The good separation be- 
tween protons and pions is clearly visible. 
Using TOF information from only one de- 
tector (e.g. hodoscope HI), the upper 
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Fig. 6. m 2 distribution for hodoscope HI up- 
per detector), 0.6 < p < 1.1 GeV/c. Log scale 
is chosen due to significant difference between 
different particles fluxes. 
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Fig. 7. m 2 distribution for hodoscope HI (up- 
per detector), 1.5 < p < 2.0 GeV/c 

bound for the separation is about 2.4 GeV/ c. 

The availability of two independent de- 
tectors HI and H2 suggests the use of 
two-dimensional distributions, as shown in 
Fig. 8. Independent constraints on the TOF 
values from both hodoscopes is the best 
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Fig. 8. Two-dimensional distribution of 
squared mass from Hodoscope (HI) and 
Preshower(H2) for the momentum region 
1.5 < P < 2.0 GeV/c. 

strategy to minimize the contamination. 
However, this is not the best way to maxi- 
mize the efficiency while maintaining a low 
contamination. In the case of two detec- 
tors a linear 'valley cut' in the plane of 
the two detector responses gives an im- 
provement [5]. It is possible to apply the 
combined constraint on the sum m 2 H1 + 
m 2 H2 . This evaluation provides an improve- 
ment in the separation between protons 
and pions - especially for higher momenta 
- and therefore extends the momentum re- 
gion for hadron identification with high ef- 
ficiency and low contamination. The iden- 
tification of kaons is limited to p < 1.5 
GeV/c, due to the very low kaon flux (see 
Fig. 6). The comparison of the pion con- 
tamination in the proton sample for indi- 
vidual constraints on the two hodoscopes 
with that for the combined constraint is 
presented in Fig. 9. 

The most important parameters for parti- 
cle identification are the efficiency for iden- 
tification of certain particle types and the 
contamination from other types. Fig. 10 
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PID TOF pions contamination in proton spectra 
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Fig. 9. Pion contamination in the proton sam- 
ple for the cases of individual constraints on 
one of the two hodoscopes with that for the 
combined constraint. 
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Fig. 10. Proton and pion identification effi- 
ciency and contamination vs momenta for the 
valley cut. 

presents both parameters for proton and 
pion identification. The valley cut allows 
identification of protons and pions up to 
2.9 GeV/c. The resulting pion sample has 
an efficiency above 98% and less than 4% 
proton contamination, while the protons 
have an efficiency of more than 85% and 
less than 6% pion contamination in the 
highest momentum bin. All analysis was 
done using semi-inclusive data samples from 
the year 1997 on polarized hydrogen tar- 
get. 



5 Conclusion 



The TOF method, as presented in this pa- 
per, allows an extension of the momentum 
range for hadron identification in HERMES 
towards lower momenta. The application 
of this method provides more statistics for 
hadrons detected within the HERMES ac- 
ceptance and extends the minimum value 
of z = -^f-, where is the hadron energy 
and v is the virtual photon energy. 
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